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ABSTRACT 

We use 16 /im, Spitzer-IRS, blue peakup photometry of 50 early-type galaxies in the 
Coma cluster to define the mid-infrared colour-magnitude relation. We compare with 
recent simple stellar population models that include the mid-infrared emission from 
the extended, dusty envelopes of evolved stars. The Ks-[16] colour in these models is 
very sensitive to the relative population of dusty Asymptotic Giant Branch (AGB) 
stars. We find that the passively evolving early-type galaxies define a sequence of 
approximately constant age ('^ 10 Gyr) with varying metallicity. Several galaxies 
that lie on the optical/near- infrared colour-magnitude relation do not lie on the mid- 
infrared relation. This illustrates the sensitivity of the Ks-[16] colour to age. The fact 
that a colour-magnitude relation is seen in the mid-infrared underlines the extremely 
passive nature of the majority (68 %) of early-type galaxies in the Coma cluster. The 
corollary of this is that 32% of the early- type galaxies in our sample are not 'passive', 
insofar as they are either significantly younger than 10 Gyr or they have had some 
rejuvenation episode within the last few Gyr. 
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1 INTRODUCTION 

Early-type galaxies (ellipticals and SOs, ETGs hereafter) are 
ancient stellar populations whose epoch of formation is re- 
lated to the large-scale structure formation in the Universe. 
This picture is supported by various studies that have used 
optical line-strength indices to determine evolutionary pa- 
rameters of ETGs in the cluster and field environments (see 
Renzini 2006 for a review). Among these studies, some au- 
thors (Bernardi et al. 2005, Clemens et al. 2006, Thomas et 
al. 2005, Sanchez-Blazquez et al. 2006a, Annibali et al. 2007) 
suggest that cluster ETGs have a luminosity weighted, mean 
stellar age 1-2 Gyr older than those in the field. Recently, 
Trager, Faber & Dressier (2008) challenged this view. An- 
alyzing the Coma cluster, they found that the 12 ETGs in 
their sample have mean single stellar population equivalent 



* This work is based on observations made with the Spitzcr Space 
Telescope, which is operated by the JPL, Caltech under a contract 
with NASA. 



ages of 5-8 Gyr, with the oldest systems being ^ 10 Gyr old. 
This average age is remarkably similar to the mean age of 
ETGs in low density environments (see e.g. Annibali et al. 
2007 and references therein). 

The emission from elliptical galaxies longward of ~ 
10 /im has a large contribution from the circumstellar dust 
around evolved stars, such as those on the AGB. This dust 
emission has been detected by Spitzer in early-type galaxies 
where it is seen as a wide emission feature near 10 /im with 
another broad feature near 18 /^m (Bressan et al., 2006). The 
spatial profile in the mid-infrared is similar to the optical 
profile (Temi et al., 2008, Athey et al. 2002), re-enforcing the 
notion that the origin is stellar rather than ISM. This cicum- 
stellar dust has also been detected directly as an extended 
envelope around nearby AGB stars (Gledhill & Yates, 2003). 
Longer wavelength infrared emission from elliptical galaxies 
is sometimes detected (Leeuw et al., 2004; Marleau et al., 
2006; Temi, Brighenti & Mathews, 2007.), but this cooler 
component is probably associated with a diffuse interstellar 
medium, rather than with evolved stars directly. 
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Figure 1. Example background subtracted, PBCD, IRS peakup images of the Coma ETGs. Contour levels are 2"/'^/100 where n=(- 
1,0,1,2. . .) mjy/pixel with a dashed contour for negative values at the lowest level. The grey-scale in each plot is scaled to the peak 
source flux. The last panel is an image of the peak-up point spread function. The cross in each panel gives the optical source position. 
The grey-scale has a square-root stretch and contour levels are 2"''^/100 where n=(l,3,5. . .) mjy/pixel. Pixels are I'.'S on a side. 



For objects older than 0.1 Gyr the mid-infrared emis- 
sion traces stars during their mass-losing AGB phase, and 
this evolutionary phase is still detected in objects with ages 
typical of globular clusters such as 47 Tucanae (Lebzelter et 
al., 2006). Because the luminosity of an AGB star depends 
on its mass and the main sequence turn-off mass decreases 
with time, the mid-infrared emission depends on the age of 
the stellar population. If all the stars in an early-type galaxy 
were created in an instantaneous burst 10 Gyr ago, the mid- 
infrared emission would be dominated by stars on the AGB 
with an initial mass of slightly less than 1 M0 (for solar 
metallicity). 

The strength of the silicate emission from the dusty en- 
velopes of evolved stars is a function of both age and metal- 
licity. As pointed out by Bressan, Granato & Silva (1998) 
however, the dependence of this emission on age and metal- 
licity is different to that in the optical. 

A colour-magnitude relation which includes a band con- 
taining the silicate emission from AGB stars therefore traces 
the AGB population as a function of the total luminosity, 
and can, in principle, be used with the optical relation to 
disentangle the effects of age and metallicity. 

Here, we use 16 images of 50 ETGs in the Coma 



cluster, made with the blue peakup detector of the IRS on 
Spitzer. We combine these with archival Spitzer IRAC im- 
ages at 4.5 /im and K-band images to construct the mid- 
infrared colour-magnitude relation. As objects at the dis- 
tance of the Coma cluster were too faint for IRS spec- 
troscopy we include 4 galaxies in the Virgo cluster for which 
we have IRS spectra. These 4 objects were selected to have 
spectra that show no evidence of activity, such as recent star 
formation or an AGN, and are intended as a 'template' for 
passive objects against which to compare the more distant 
objects of Coma. 

The Coma cluster is both very rich and dynamically 
relaxed, and contains some of the the most massive galaxies 
in the local Universe. As such, it is an ideal place to study 
the star formation history of ETGs. 



2 DATA ANALYSIS AND RESULTS 
2.1 Mid-infrared photometry 

The blue peakup images (see Fig [T] for examples) are back- 
ground dominated, with source fluxes only a few percent 
above the background level. In addition, the background 
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Figure 2. Relation between the effective H-band radii (1.65 /im, 
Gavazzi et al. 2000) and the 16 /xm effective radii. The latter have 
been estimated by convolving a model r^/* de Vaucouleur's law 
with the IRS blue peakup PSF (see text). Crosses are ellipticals 
and diamonds are lenticulars (classified either as SO/E, SO, SBO, 
or so/a). The solid line is a least squares fit to the data whereas 
the dotted line corresponds to equal Re at both wavelengths. The 
horizontal dot-dashed line corresponds to the radius containing 
half of the encircled energy for the IRS blue peakup PSF. 

shows fluctuations on spatial scales of several pixels. This 
presents a serious problem in the consistent determination 
of integrated fluxes for sources of different angular sizes. 
The standard technique of modeling the background as a 
smoothly varying light distribution, and then subtracting 
it from the image, tends to over-correct in the vicinity of 
larger sources, simply because the source size is similar to 
that of the background fluctuations. In order to avoid this 
problem we have opted to measure colours, rather than in- 
tegrated fluxes, and to do so within a fixed angular radius 
of 12". The flux within this radius is background corrected 
by subtracting the median pixel value within an annulus 
just outside this aperture. The error on the derived flux was 
taken as the quadrature sum of the rms variation within the 
annulus and the calibration error, which was taken as 5% 
for the blue peakup data. 

Before applying the same procedure to IRAC channel ^ 
(4.5 /im) and 2MASS Ks-band (2.2 /im) images, these were 
convolved to the resolution of the IRS-peakup data, in order 
to remove colour gradients due to the differing instrumental 
resolutions. The measured aperture fluxes were then used to 
derive the colours for each object; these colours are presented 
in Table [U 

Aperture photometry for the 4 Virgo elliptical galaxies 
was carried out in an analogous manner, but rather than use 
a circular aperture of radius 12", elliptical apertures of semi- 
major axis 72" were used. This difference takes account of 
the fact that the Virgo cluster is 6 times closer than Coma 
cluster and that the images are well-resolved. There is a 

^ Calibration errors taken as 10% , although the real calibration 
errors may be smaller. 



possibility that the use of a larger aperture in the presence 
of a fixed psf will introduce a systematic error in the flux 
determinations between the Coma and Virgo galaxies. In 
order to verify that this was not a siginficant source of error 
we constructed a number of pairs of simulated galaxy images 
with de Vaucouleur profiles, differing in their effective radii 
by a factor of 6. We then convolved each with the IRS blue 
peakup psf and measured the flux within apertures of 12" or 
72" accordingly. The differences in the retrieved fluxes did 
not exceed 10%, and were typically less than the photometric 
errors. 

All aperture photometry was performed using a custom 
pipeline written in idl and use was made of the IDL As- 
tronomy User's Library (Landsman 1995). 

2.2 Estimating the effective radii 

A direct estimate of the effective radius based on an im- 
age with a high background level and fairly low signal-to- 
noise ratio would be rather sensitive to any error in the 
subtraction of the background level. In order to minimize 
this problem, rather than finding the radius containing half 
of the integrated flux, we adopted the following procedure. 
We flrst constructed a model galaxy with a de Vaucouleur's 
r^'^* law and convolved it with the IRS blue peakup point 
spread function (PSF). We either assumed that the intrinsic 
ellipticity of each galaxy at 16 /im was the same as the op- 
tical value as given in "Goldmine'13 or was circular, where 
no value was available. After convolution with the PSF we 
measured the radial profile in circular apertures and com- 
pared this to the same apertures measured on the image. 
The derived value of the effective radius, i?e, is that of the 
model galaxy whose PSF-convolved FWHM was equal to 
that of the image profile. These values are given in Table [1] 
along with effective radii measured at 1.65 /.im by Gavazzi 
et al. (2000). The error on the estimate of the effective radii 
was estimated using a Monte-Carlo approach. To the model 
galaxy, convolved with the IRS blue peakup PSF, Gaus- 
sian noise was added at the level present in the real images. 
For 200 realizations of this noise, the FWHM was measured 
using the same procedure as for the real images. The distri- 
bution of these values (as a function of S/N ratio) was used 
as an estimate of the error on the derived values of Re. 

A comparison between the H-band and 16 /xm effective 
radii is shown in fig. [5] Clearly, most galaxies in Coma are 
barely resolved by Spitzer-IRS, however, for those that are 
clearly resolved their effective radii are similar to those in 
the near-infrared. There is no evidence, for example, that 
the 16 /im emission is derived from an unresolved central 
source. Similar effective radii in the near and mid-infrared 
are expected if the mid-infrared emission is stellar in origin, 
such as found by Temi et al. (2008). 

2.3 The mid-infrared colour-magnitude diagram 

In Figure 13] we show the mid-infrared colour-magnitude di- 
agram of the Coma cluster. The Ks-[16] colour is shown as 
measured within a fixed aperture of 12 arcsec. 

^ |http : //goldmine ■ mlb ■ Inf n . It/ 1 , Gavazzi et al. (2003) 
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Table 1. The integrated 16 jxra fluxes and sizes of the sample ETGs. The flux errors are the quadrature sum of a term due to the 
uncertainty in deriving the background sky level and a 5% cahbration uncertainty. The effective radii have been calculated by convolving 
a model r^/* de Vaucouleur's law with the IRS blue peakup PSF (see text). 



The distribution of Kg- [16] colour shows a steady red- 
dening from the faintest magnitudes observed (Ks ~ —22) 
to Ks ^ —24. In this magnitude range the dispersion in 

Kg — [16] is close to a magnitude with values as high as 2.2 
near Ks ~ —24. At brighter magnitudes the colours show 
little variation and low dispersion with Ks — [16] ~ 1.3. 

We also plot the position of four selected passive ETGs 



in the Virgo cluster, selected from the sample of Bressan et 
al. (2006). These galaxies have nuclear MIR spectra whose 
main characteristic is the presence of the silicate bump above 
10 /im, indicating a passively evolving stellar population. 
These galaxies have also been observed in the IRS Peakup 
mode to secure a link between photometric and spectro- 
scopic observations. 
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Figure 3. Mid-infrared colour-magnitude relation of the Coma 
cluster for fixed aperture Ks-[16] colours. Objects are identified 
according to their designations in Dressier (1980). The 4 ETGs 
in the Virgo cluster are shown with coloured symbols. For each 
of the Virgo galaxies the mid-infrared colours within the central 
10" radius, as well as that of the whole galaxy excluding this 
central aperture are shown, as upward and downward pointing 
symbols, respectively. The assumed distance moduli to the Coma 
and Virgo clusters are 35.1 and 31.3 respectively. 
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Figure 4. As fig. [3] but for photometry with circular radii equal 
to 2iJe where Rc is the effective radius of the galaxy measured 
in the H-band (Gavazzi et al., 2000). For the 4 Virgo galaxies a 
radius of only Rc was used due to the limited size of the pcakup 
images. 



The presence of stellar population gradients within 
ETGs (Annibali et al. 2007), together with the existence of 
a general increase of the effective radius with the total mag- 
nitude of the galaxy, could, however, introduce a spurious 
effect, since a fixed aperture measurement would probe stel- 
lar populations in relatively more central regions for brighter 
galaxies. 

To quantify this effect we repeated the measurement 
within a fixed relative aperture of radius ~2 Re, where i?e 
is the effective radius of the galaxy. The results are shown 
in fig. m though the lack of measured H-band effective radii 
introduces a cut at magnitudes fainter than Ks ~ —23. In 
this figure the flattening at magnitudes brighter than Ks ~ 
—24 becomes more evident and actually the slope of the 
colour-magnitude relation is reversed for the brightest Coma 
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Figure 5. The optical/NIR colour-magnitude relation of the 
Coma cluster using the photometry of Bower, Lucey &; Ellis 
(1992). Coma ETGs are shown as squares whereas the passive 
ETGs in the Virgo cluster are shown as triangles. Galaxies with 
a Ks-[16] excess are indicated with a cross. 



galaxies, while there are only small differences at magnitudes 
fainter than Ks ~ —24. 

The position of the four passive Virgo galaxies in fig. [3] 
and that of the Coma galaxies brighter than Kg ^ —24 define 
a tight relation that can be extended down to the faintest 
magnitudes. By analogy with optical/NIR colour-magnitude 
relations and because of the presence of the four Virgo ETGs 
we argue that this line represents the MIR colour-magnitude 
relation of passively evolving, ETGs in the Coma cluster. 
ETGs significantly above this line have an "excess" of 16 /im 
emission with respect to the passive sequence. 

The presence of a population of galaxies which appear 
anomalously red in the mid-infrared colour-magnitude di- 
agram is not mirrored in the optical/near-infrared, fig. (5] 
The V-K colour shows a steady increase towards brighter 
magnitudes with a slight fiattening for objects brighter than 
K ~ —24. This, in fact, is very similar to what we identify 
above as the passive colour magnitude relation in the mid- 
infrared. The implications of this difference are discussed 
further below. 



3 DISCUSSION 

From an extrapolation of the NIR stellar photospheric con- 
tinuum into the mid-infrared Bressan et al. (2006) have es- 
timated that the 16 ^m emission in passive, ETGs has ap- 
proximately equal contributions from stellar photospheres 
and dust emission. The most likely origin of the latter emis- 
sion is from dusty circumstellar envelopes around evolved 
stars, such as the mass losing Asymptotic Giant Branch 
(AGB) stars. The exact proportions depend on parameters 
such as the age and metallicity of the stellar population. 
Because the photospheric contribution shows little varia- 
tion for a fixed K-band magnitude, the mid-infrared colour- 
magnitude relation should depend strongly on the popula- 
tion of dusty AGB stars as a function of galaxy luminosity. 
In this scenario, those ETGs that show an excess of 16 /im 
emission have the largest luminosity-weighted dusty AGB 
population. The evolution of this stellar population should 
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therefore be critical to the interpretation of the mid-infrared 
colour-magnitude diagram. 

Another possible origin of an excess of emission at 
16 fim could be the contribution of a central AGN, as seen in 
the spectrum of NGC 4486 (Bressan et al. 2006). However, 
there is no evidence of broad emission lines in the avail- 
able SDSS spectra of those galaxies which lie well above the 
'passive' colour-magnitude relation in fig. |3l There is also no 
evidence that these galaxies have smaller efi'ective radii at 
16 /im (as shown by NGC 4486, that appears unresolved in 
the SL and LL IRS apertures). 

3.1 Influence of the inter-stellar environment 

Before considering the colour-magnitude relation in terms 
of stellar evolution, we first consider the infiuences that the 
inter-stellar environment might have on the mid-infrared 
colours. In nearby AGB stars a significant fraction of the 
16 /im fiux originates from extended, dusty, circumstellar 
envelopes. It is thus important to understand how robust 
these envelopes are to different external perturbations, such 
as the ram-pressure of the interstellar medium (ISM) in the 
central regions of massive ellipticals. 

The circumstellar envelope of an AGB star expands at 
a typical velocity of ~ 10 kms~^ which, at the condensation 
radius, Rc ~10^''cm (Bressan et al. 1998, Granato & Danese 
1994) is of the order of the escape velocity. The MIR emis- 
sion, f(!^), in AGB stars, has a broad maximum between 10 
and 20 /im suggesting a dust temperature Tdust ~ 300 K. 
Inspection of our dusty circumstellar models shows that dust 
reaches this temperature only in a thin internal region of the 
envelope, whose size is a few tens of Rc- At 10 kms~^ the 
crossing time of this region is only ICQ yr. Thus, the MIR 
emission from AGB stars is a short lived phenomenon: after 
about 100 yr, dust is already so cool that its contribution 
in the MIR spectral region becomes negligible. However, in 
the central regions of an ETG, the star moves at the dis- 
persion velocity through the tenuous ISM, and we need to 
determine whether or not the circumstellar envelope is able 
to maintain its structure in the inner few tens of Re- 
in order to estimate the radius at which the organized 
structure of the envelope is destroyed by its interaction with 
the ISM, we equate the kinetic energy density in the circum- 
stellar wind (~ l/2pcv'^) to the kinetic energy density of 
the ISM, due to the stellar dispersion velocity, (~ l/2pi*ag). 
Using a stationary wind the circumstellar gas density is 



PC 



M 



4:TVr'^Va 



(1) 



while the ISM density is evaluated from the mass that has 
been already lost by stars 



Pi = 



dtdM, 



X t X p. 



X t X pt 



(2) 



In the second equation d? M /dtdMt, ~10~^^yr~^ is the typ- 
ical gas deposition rate per unit time and per unit stellar 
mass from the old stellar population. This is obtained from 
integration of the mass loss rate along our isochrones of 
old stellar populations (see also Tinsley 1973). The average 
stellar density, p*, is evaluated from deprojection of the De 
Vaucouleurs law (Young 1976) and depends on the assumed 
position within the galaxy. 



Applying the condition of virial equilibrium and the ob- 
served relation between the total galaxy mass and the ve- 
locity dispersion (Cappellari et al. 2006) we find 



r/r-c 



200i 



POltCy 



(3) 



In this equation, Tc is the condensation radius, ~ 1.66 x 
lO^'^ y^L/LQ , for a silicate dust mixture (for graphite ic is 
about half this value and t/tc must be doubled); AIs is the 
stellar mass loss rate in lO~*M0/yr; vio is the wind velocity 
in units of 10 kms~^; L3 is the luminosity in units of IO^Lq; 
as is the stellar velocity dispersion in units of 300 kms"'^; Ve 
is the effective radius of the galaxy; poi is the dimensionless 
average stellar density with respect to its value ai r/ve — 0.1 
(poi =10.5, Young, 1976); M12 is the galaxy mass in units 
of IO^'^Mq and tcyr is the time in Gyr. 

Thus, only for very low mass-loss rates (lO^^^MQ/yr) or 
in the very central regions of a galaxy - p(r jr^ — 0.01) — 357 
- may the structure of the wind be destroyed at a few con- 
densation radii. We may thus conclude that, in all relevant 
cases, the circumstellar wind is not destroyed by its interac- 
tion with the environment even in the cores of ETGs. 

Once the circumstellar envelope mixes with the ISM, 
dust grains may be destroyed via sputtering by the hot (~ 
10^ K) ISM gas. This process, however, will not affect the 
MIR emission, because, as discussed above, when the dust is 
released to the ISM it is already too cool to emit in the MIR 
spectral region. Furthermore, the timescale for dust grain 
destruction via sputtering, even in the centre of the Coma 
cluster, is ~ 10* yr (Dwek, Rephaeli & Mather, 1990), much 
larger than the lifetime of the newly formed hot dust within 
the circumstellar envelope (~100 yr). 

In summary, we conclude that the MIR emission from 
dust grains in AGB stars is a short lived phenomenon 
(~ 100 yr) happening in the inner, denser region of the cir- 
cumstellar envelopes. This emission cannot easily be modi- 
fied by interaction with the environment, so that the envi- 
ronment does not have a direct infiuence on the mid-infrared 
colour-magnitude relation. 



3.2 Stellar Evolution 

As it appears that the form of the MIR colour-magnitude 
relation is not influenced directly by some particular envi- 
ronmental effect, it must be explained via the properties of 
the underlying stellar populations. 

In Figure |S] we compare our observations with the pre- 
dictions for the integrated colours of simple stellar popula- 
tions (SSPs). 

As discussed in Bressan et al. (1998) both the integrated 
optical and MIR colours of SSPs suffer from the degeneracy 
between age and metallicity, but in the opposite sense. While 
optical features (colours or narrow band indices) of the same 
strength can be obtained by an anti-correlated variation of 
age and metallicity, producing the typical banana-shaped 
solutions in the age-metallicity diagrams (Annibali et al. 
2007), an equal intensity of the 10pm bump may be ob- 
tained by a correlated variation of the age and metallicity. 

For this reason, the effects of age and metallicity varia- 
tions are best separated in a mixed optical-NIR-MIR colour- 
colour plot, such as that in Fig. [S] 
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Figure 6. Optical NIR MIR colour-colour diagram. The horizon- 
tal axis is the V-K colour from the precision photometry of Bower, 
Lucey & Ellis (1992) while the vertical axis is the Ks-[16] colour 
measured within fixed apertures of 12". Squares are Coma ETGs 
and triangles are the four selected passive ETGs in the Virgo clus- 
ter. Galaxies with a Kg- [16] excess are indicated with a cross. SSP 
models for the indicated ages and metallicity are computed fol- 
lowing the scheme of Bressan et al. (1998). Models without dusty 
AGB envelopes collapse to the almost horizontal dashed lines at 
the bottom of the diagram. The dotted line represents the locus of 
the Kurucz model atmosphere with low metallicity (the effect of 
metallicity is very small however) and the asterisk is a COMARCS 
low temperature atmosphere model as calculated by Aringer et 
al. (2008, in prep.) The position of the globular cluster, 47-Tuc, 
is indicated by the solid hexagon in the lower left. 
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Figure 7. Optical SDSS spectra of selected Coma cluster galax- 
ies. Galaxies with a Kg- [16] excess are shown in green and are 
indicated with an asterisk. 
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Figure 8. Ks-[16] colour as a function of the projected distance 
from the cluster centre, defined as the CD galaxy D148 (asterisk 
on y-axis). The bottom panel shows the cumulative distribution 
for ETGs with Ks-[16] > 1.5 and Ks-[16] < 1.5. The number of 
objects in each sub-set is indicated in parentheses 



In Fig. [6] the horizontal axis is the V-K colour ob- 
tained from the precision photometry of Bower, Lucey & 
Ellis (1992), while the vertical axis is our measured Kg- [16] 
colour. The K-band magnitudes on the two axes therefore 
refer to slightly different central filter wavelengths and aper- 
tureslfl Though the V-K colours of Bower, Lucey & Ellis are 
provided only for a sub-sample of our galaxies (30 out of 50 
objects), we prefer to use their values to preserve the homo- 
geneity of the data. Squares are the Coma ETGs while the 
triangles refer to the four Virgo galaxies. 

New SSP models, for the indicated age and metallic- 
ity, have been computed following the scheme of Bressan et 
al. (1998). The only difference with respect to these previ- 
ous models is a minor revision of the mass loss rate during 
the AGB phase, introduced by the calibration with new V- 
K integrated colours of populous star clusters in the Large 
Magellanic Cloud, and the adoption of the new empirical 
library of stellar spectra (Sanchez-Blazquez et al., 2006b). 

A few caveats concerning models should be kept in mind 
because they are relevant to the present investigation. 



^ The photometry of Bower, Lucey & Ellis was measured in an 
aperture of radius 8'.' 5 and corrected to a radius of 5'.' 5. The aper- 
ture correction applied was 0.03 mag for V-K. Assuming a cor- 
rection of similar magnitude would be required to correct to our 
radius of 12", the V-K values should shift by -0.03 mag. This 
would not alter any of our interpretation. 



8 Clemens et al. 



First of all, it is worth recalling that in Bressan et al. 
(1998) the effects of possible circumstellar dust in the red gi- 
ant branch (RGB) phase (i.e. before central helium burning 
in the horizontal branch in the case of old populations) have 
not been considered. Recent claims, that RGB stars also 
show a significant MIR excess not dependent on their lumi- 
nosity (Origiia et al. 2007), should be regarded with great 
care because of possible strong crowding effects (Boyer et al. 
2008). 

Furthermore, old stellar populations in the galaxies we 
are considering could reach super-solar metallicity and, as 
such, may not be well represented by those in our own 
galaxy. Very little is actually known on the efficiency of 
the mass-loss rate at high metallicity (Carraro et al. 1996). 
However, van Loon et al. (2008) do not find evidence of a 
MIR excess in the RGB stars of NGC 6791, a super metal- 
rich globular cluster. Also, NGC 4649 and NGC 1399, the 
two ETGs with the largest UV upturn and high Mg2 index 
(Bertola et al. 1995), that have been interpreted as the sig- 
nature of very metal-rich stars (Bressan, Chiosi & Fagotto 
1994), have Spitzer IRS spectra similar to those observed 
in our passive Virgo ETGs (Bressan et al. in preparation). 
Thus, there seems to be no evidence that metallicity has a 
direct effect on mass loss. 

As already anticipated, variations of age and metallicity 
are very well separated in this colour-colour diagram, espe- 
cially at high metallicity, with iso-metallicity lines being al- 
most horizontal and coeval SSPs moving on almost vertical 
lines. In contrast, models without dusty AGB envelopes col- 
lapse to the almost horizontal dashed lines at the bottom of 
the diagram, independent of their ages. The effects of dust 
disappear at ages larger than 13 Gyr, where the AGB phase 
is practically absent. 

We have extended the photospheric stellar spectra of 
MILES in the NIR/MIR region by means of matched 
NEXTGEN models (Hauschildt et al. 1999). In addition, 
the effects of dusty envelopes have been included following 
the prescriptions described in Bressan et al. (1998). Thus, 
the bottom horizontal dashed lines (models without dusty 
AGB envelopes) superimpose on single NEXTGEN spectra 
at varying effective temperature. The effect of varying the 
metallicity is almost negligible in the Kg- [16] colour. To illus- 
trate the difference introduced by adopting different stellar 
spectra we have marked with an asterisk the position of a 
COMARCS low temperature atmosphere model as calculated 
by Aringer et al. (2008, in prep.). Finally, the hexagon to- 
wards the bottom left of the diagram (age ~ 12 Gyr and Z 
^ 0.004), marks the position of the globular cluster 47 Tuc 
whose semi-empirical integrated colours have been derived 
by Bressan et al. (2006b). 

Fig. [6] shows that galaxies that follow the MIR colour- 
magnitude relation are almost coeval, with an age of about 
10 Gyr and a metallicity that grows by a factor of about 
two. It is interesting to note that the four Virgo galaxies 
are as old as the oldest Coma galaxies. None of the galaxies 
appear to be as old as the globular cluster 47 Tuc. 

Galaxies with a Kg- [16] excess, cluster around the lines 
of 8 Gyr and 6 Gyr thus appearing younger by about 2-4 
Gyr. 

In order to clarify the different sensitivity to stellar pop- 
ulations of optical and MIR colours, we have reproduced, in 
Fig. [5l the optical/NIR colour-magnitude relation for the 



galaxies shown in Fig.|6] Galaxies with a Kg- [16] excess (i.e. 
clustering around and above 8 Gyr in Fig. [Sjl are plotted 
with a cross. Fig. [S] shows that in spite of a well detected 
Ks-[16] excess, the latter galaxies lie on the optical/NIR 
colour-magnitude relation. Moreover, a significant fraction 
of these galaxies are classified as ellipticals. 

3.3 Indirect influence of cluster environment on 
the Ks-[16] colour ? 

The fraction of galaxies showing a Kg- [16] colour excess with 
respect to the oldest systems, in Fig. [6l is about 48%. If we 
consider the total sample (Fig. [3} and add four more galaxies 
with a clear Ks-[16] excess with respect to the CM relation, 
the fraction is 32 ± 8%. 

A possible interpretation, discussed in the previous sec- 
tion, is that these systems have average younger ages, since 
a younger stellar population is expected to have a redder 
MIR emission by virtue of the larger fraction of dusty AGB 
stars. However, galaxies are complex systems that may not 
be well represented by simple stellar populations and cau- 
tion must be used in the interpretation of the Kg- [16] colour 
excess. It does not necessarily point to a younger global age 
of a passively evolving system. 

For example, in the case of the Virgo ETC, NGC 4435, 
Panuzzo et al (2007) have shown that the MIR emission is 
dominated by a nuclear starburst which is practically invis- 
ible in the optical and NIR. Its IRS low resolution spectra 
are dominated by PAH, atomic and molecular emission fea- 
tures that are typical of star forming galaxies. At 16/im the 
emission is not resolved, indicating that the emission arises 
from the nucleus; which indeed appears as a dusty disk in 
HST images. In this rejuvenation episode (about 200 Myr 
ago) less than 1% of the total galaxy mass was formed. From 
Panuzzo et al. (2007) we estimate, for NGC 4435, a Ks-[16] 
colour ~2.8 mag in the central 5 arcsec radius aperture, 
and a total Ks-[16] colour ~2.1 mag and total V-K ~3.1. In 
Fig. |6] the latter values would put NGC 4435 in the region 
around 6 Myr and Z=0.02. Indeed, NGC 4435 was selected 
for IRS observations as one of the galaxies that make up the 
Virgo colour-magnitude relation. From the analysis of this 
sample (18 galaxies) Bressan et al. (2006) concluded that 
the fraction of early-type galaxies that show signatures of 
recent activity is 17%. This figure is significantly lower than 
the fraction of Coma ETGs with a Kg- [16] colour excess. 

We have also extracted, from data release 6 of the Sloan 
Digital Sky Survey (SDSS) the available spectra (16 out 
of 50) of our Coma galaxy sample. In Fig. [7] we show the 
3700A-4200A spectral region, around the Call lines. Galax- 
ies are sorted by decreasing brightness and green curves rep- 
resent objects with a Ka-[16] colour excess. In three objects 
(D103, D118 and D155) with a Ks-[16] excess, the Call H 
line is as deep as the K line, indicating contamination by A 
type stars, i.e. a recent rejuvenation episode. In the other 
galaxies the Call doublet appears normal. We recall, that 
for solar metallicity, an inversion of the Call doublet starts 
appearing only at ages younger than 2 Gyr (Longhetti et 
al., 1999). 

Although we see in Fig. |8] that the Ks-[16] colour is un- 
correlated with the position in the cluster we note that if we 
identify those galaxies with an 'excess' of 16 urn emission 
as those with Kg — [16] > 1.5 (see fig. [3)l then the mean 
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projected distance from the cluster centre of this popula- 
tion is smaller than for those galaxies with Ks — [16] < 1.5. 
The mean projected distance from the cluster centre for the 
16 /xm excess galaxies is 0.17° whereas that for the oth- 
ers is 0.32° (the KS-probability that the 2 distributions are 
the same is 1.6%). Because it seems that the cluster envi- 
ronment does not directly influence the mid-infrared colour- 
magnitude diagram this suggests that the star formation his- 
tory has been modified by the environment; which in turn 
influences the galaxy colours. The sense of this influence is 
that galaxies with redder Ks-[16] colour (younger luminosity 
weighted stellar ages, flg.O are preferentially found at small 
cluster-centric radii. 

Miller, Neal & Owen (2002) have investigated the dis- 
tribution of both AGN and star-forming galaxies as a func- 
tion of cluster-centric distance. They find that dusty star- 
forming galaxies have a more centrally concentrated distri- 
bution than normal star-forming galaxies, with AGN being 
more centrally concentrated than both. Though there is no 
evidence of emission lines in the SDSS spectra, we cannot ex- 
clude that the MIR excess is due to an obscured nuclear star- 
burst like that found in NGC 4435. A possible way to check 
this possibility is by means of existing IRAC 4.5/im and 8/im 
observations (sampling the strength of the 7.7 PAH feature) 
and MIPS 24/xm (sampling the presence of hot dust). This 
work is underway. 



3.4 Implications 

The MIR colour-magnitude relation indicates that the old- 
est Coma early-type galaxies define a sequence of increasing 
metallicity at constant age [Z ~ 0.01 — 0.03, i ~ 10 Gyr). 

While this conclusion has already been suggested by the 
analysis of the optical/NIR colour-magnitude relation, it is 
only with the use of the complementary MIR observations 
that the alternative solution (invoking a synchronization of 
the epoch of formation with total brightness, see Bower, 
Lucey & Ellis, 1992) can be excluded. 

Indeed, our analysis of the mixed optical-NIR-MIR 
colour-colour diagram clearly shows that some of the galax- 
ies that populate the optical colour-magnitude relation may 
have average ages that are even 40% younger than the old se- 
quence. In those galaxies, however, a larger metallicity may 
redden the colours and compensate for the younger age in 
an optical diagram. This illustrates the importance of the 
MIR colour in breaking the age-metallicity degeneracy. 

Our conclusion is in very good agreement with the anal- 
ysis of narrow band indices of ~ 4000 ETGs carefully se- 
lected from the SDSS survey (Clemens et al. 2006) and re- 
cently extended to ~ 14000 objects with the DR6 catalogue 
(Clemens et al. 2008 in prep.). These authors find that the 
luminosity-weighted ages of the whole sample start young 
at a low velocity dispersion, rise until a ~ 2.32 and fiatten 
at larger a. In contrast, the metallicity shows a continuous 
increase with increasing o. 

Our result, that a large fraction of ETGs in Coma pop- 
ulate a genuine old and coeval sequence, contrasts somewhat 
with the recent study of ETGs in the Coma cluster by Trager 
et al. (2008). Based on optical spectra these authors also 
find that, in spite of a wide range in mass, 10/12 objects are 
consistent with the hypothesis of a coeval population. How- 



ever, the age they find for this population, 5.2 ± 0.2 Gyr, is 
significantly younger than our value. 

We notice that this discrepancy is unlikely to be caused 
by a significant offset between ages derived with the two 
methods. Indeed, there are four objects in common, for 
which Trager et al. find the following ages: D129 (GMP3329) 
7.9 Gyr, D157 (GMP3484) 7.8 Gyr, D155 (GMP3367) 4.5 
Gyr and D133 (GMP3639) 3 Gyr. 

D129 and D157, the two old objects, faU on our MIR 
colour-magnitude relation, for which we derive only sligthly 
older ages. On the other hand D155 and D133, the two young 
objects, lie well above this relation and, given their Ka-[16] 
colour we would conclude that they are possibly affected by 
some rejuvenation episode. 

Considering that, contrary to our analysis, optical stud- 
ies are affected by the age-metallicity degeneracy causing 
(large) uncertainties that are difficult to quantify, we be- 
lieve that there is good agreement for the few objects in 
common. 

It is more difficult to explain why we find that 68% of 
the galaxies in our sample populate the old sequence, while 
Trager et al. find that only 25% are old. If their sample is 
not biased toward young objects and given the sensitivity 
of the Ks-[16] colour to rejuvenation events, we suspect that 
the age-metallicity degeneracy is still to blame. 

Bregman et al. (2006) find a very similar discrepancy 
between the age of ETGs determined by optical line indices 
and mid-infrared Spitzer spectra that include the silicate 
emission feature near 10 fim. They also find mean ages of 
~ 10 Gyr using mid-infrared data. 

It is also interesting to note that the upper 2.5 mag of 
the MIR colour-magnitude relation is quite narrow. On the 
contrary, in the lower 2.5 mag, below Mk ~ —24.5, there 
is a significant dispersion. This shows that the rejuvena- 
tion events (of any kind) have not affected the most massive 
ETGs. 



4 CONCLUSIONS 

We present 16 pim, Spitzer-IRS, blue peakup images of a 
sample of 50 ETGs in the Coma cluster. We compare these 
with archival IRAC images at 4.5 //m and 2MASS, Ks band 
images at 2.2 /im. 

We make the following conclusions. 

• Our IRS blue peakup images show no evidence that the 
16 /im emission is anything other than stellar in origin 

• The region within a dusty AGB star envelope where 
dust is hot enough to emit strongly at 16 is not vulner- 
able to environmental effects such as the ram-pressure of an 
ISM wind or dust sputtering by hot gas. Dust grains only 
spend ^ 100 yr near ^ 300 K and this region is sufficiently 
dense to survive against ram-pressure disruption. 

• We identify the mid-infrared colour-magnitude relation 
of passively evolving ETGs as the lower envelope of the 
galaxy distribution in the Ks - [16] vs Ks plane; that is, the 
minimum value of Ks - [16] at a given Ks magnitude. ~68% 
of the galaxies in our sample lie on this colour-magnitude 
relation. These galaxies cannot have had any episode of star 
formation accounting for more than ^ 1% of the total stellar 
mass within the last few Gyr. These are genuinely 'passively 
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evolving' objects. The remaining objects are either signifi- 
cantly younger than 10 Gyr or have undergone a rejuvena- 
tion event in the recent past. This result is at odds with the 
most recent estimate of the fraction of old objects based on 
optical spectroscopy (Trager et al. 2008). 

• We construct the mixed optical-NIR-MIR two colour di- 
agram and, by means of updated simple stellar population 
models, we show that the addition of mid-infrared data al- 
lows a much better separation of the effects of age and metal- 
licity, which are rather degenerate in either the optical or 
mid-infrared when taken in isolation. In this plane galaxies 
populating the colour-magnitude relation trace a sequence of 
varying metalUcity at approximately constant age. Although 
this conclusion was already consistent with the optical-NIR 
colour-magnitude relation, a correlation between age and 
metallicity could equally well explain the relation. Indeed, 
comparison with the optical-NIR colours shows that a num- 
ber of galaxies that lie on the optical-NIR relation are signif- 
icantly displaced from the mid-infrared relation, with redder 
Ks - [16] colours. The mid-infrared colour-magnitude dia- 
gram therefore shows a sequence of metallicity for old, pas- 
sive galaxies, with younger objects displaced towards redder 
Ks - [16] colours. 

• The oldest elliptical galaxies in our sample have lumi- 
nosity weighted, mean stellar ages of 10.5 Gyr and metallic- 
ities within a factor of two of the solar value. No galaxy in 
our sample is as old or metal poor as the globular cluster 47 
Tuc. 

• Although the addition of the Ks - [16] colour allows 
us to identify objects with significantly younger, luminosity 
weighted, mean stellar ages, we cannot distinguish between 
genuinely 'young' objects and those that have undergone a 
minor rejuvenation event. However, given that even a period 
of recent (last few Gyr) star formation that accounts for 
less than 1% of the total stellar mass will shift a galaxy off 
the mid-infrared colour-magnitude relation, the latter option 
seems far more likely. Unambiguous resolution of this issue 
will require the infrared spectroscopic capabilities of future 
space observatories. 

• There is evidence that those galaxies that have an ex- 
cess in the Ks - [16] colour are found preferentially at smaller 
cluster-centric radii. As the interaction between the dusty 
AGB star envelopes and the ISM does not directly effect 
the 16 fim emission, the excess may be caused by "rejuve- 
nation" episodes induced by the cluster environment. 
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